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ABSTRACT 

We search for the maximum oxygen abundance in spiral galaxies. Because this maxi- 
mum value is expected to occur in the centers of the most luminous galaxies, we have 
constructed the luminosity - central metallicity diagram for spiral galaxies, based on 
a large compilation of existing data on oxygen abundances of Hn regions in spiral 
galaxies. We found that this diagram shows a plateau at high luminosities (-22.3 < 
Mb < -20.3), with a constant maximum value of the gas-phase oxygen abundance 
12+log(0/H) ~ 8.87. This provides strong evidence that the oxygen abundance in 
the centers of the most luminous metal-rich galaxies reaches the maximum attainable 
value of oxygen abundance. Since some fraction of the oxygen (about 0.08 dex) is 
expected to be locked into dust grains, the maximum value of the true gas+dust oxy- 
gen abundance in spiral galaxies is 12+log(0/H) ~ 8.95. This value is a factor of ~ 2 
higher than the recently estimated solar value. Based on the derived maximum oxygen 
abundance in galaxies, we found the oxygen yield to be about 0.0035, depending on 
the fraction of oxygen incorporated into dust grains. 
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1 INTRODUCTION 

It is well known that there is a large scatter in the heavy 
element content of galaxies. This is due to the fact that dif- 
ferent galaxies evolve chemically at different rates. Galactic 
winds with different efficiencies can also make a significant 
contribution to the scatter in metallicity in low-mass dwarf 
irregular galaxies. The oxygen abundance in the interstel- 
lar gas is usually used as a tracer of metallicity in late-type 
(spiral and irregular) galaxies. A question of great interest 
in the chemical evolution of galaxies is that of the maxi- 
mum value of the observed oxygen abundance. Is there such 
a value? And is this maximum value equal to the maximum 
attainable value of the oxygen abundance in galaxies? 

The latter is defined by the stellar oxygen yield, i.e. the 
mass of oxygen freshly synthesized and ejected by a gen- 
eration of stars relative to the mass locked up in low-mass 
stars and compact remnants. It can be derived with in the 
fram ework of chemical evolution models of galaxies |Pagell 
1 19971 ). The closed-box model which neglects mass exchange 
between a galaxy and its environments gives the maximal 
upper bound to the metallici ty of the gas for a given gas 
mass fraction (|Edmundslll990l ). In practice, the procedure is 



often inverted, i.e. the measured chemical compositions of 
galax ies are used to estimate empirical ly the oxygen yield 
(e.g. lGarnettJ[200l : iPilvugin et al.l l2004). This is so because 
the theoretical value of the oxygen yield is not well known, 
due to uncertainties in both the oxygen production by stars 
of different masses and metallicities and the parameters of 
the initial mass function of stars, the relative birthrates of 
stars with different initial masses. Thus the uncertainty in 
the stellar oxygen yield prevents an accurate determination 
of the maximum attainable value of the oxygen abundance 
in galaxies through chemical evolution models of galaxies. 

Can we determine that maximum value by observa- 
tions? Oxygen abundance s in the most metal-rich spira l 
galaxies in the samp le s of IVila-Costas fc Edmunds! (1 199*3) 



IZaritskv et ail (| 19941 ): iGarnett et all l| 19971) 



, , Ivan Zee et al 

ljl998t)) have been estimated previously by IPilvugin et al 



l|2006al ). Those authors found the maximum gas-phase oxy- 
gen abundance in the central regions of those spiral galaxies 
to be 12+log(0/H) ~ 8.75, suggesting that this value may 
be an upper limit to the oxygen abundances in spiral galax- 
ies. 

Here we carry out a search for the maximum oxy- 
gen abundance using a considerably larger sample of spi- 
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ral galaxies. The strategy of ou r search is base d on th e two 
following considerations. First, iLequeux et al.l (|l979l ) have 
found that, for irregular galaxies, the oxygen abundance 
correlates with the total galaxy mass, in the sense that the 
higher the total mass, the higher the heavy element content. 
Since the galaxy mass is a poorly known quantity, the lumi- 
nosity - metallicity (L - Z) relati on is often used ins t ead o f 
the mass - metallicity relation. iGarnett fc Shields! (|l987l ) 
have indeed found that spiral disk abundances correlate 
well with galaxy luminosities. This luminosity-metallicity 
correlation has been confirmed in many studies of spi- 



ral g a laxies (|Vila-Costas fc Edmunds! [1993 : Zaritsky ct al 



1994; iGarnetj 20021; IPilvugin et al.l |2004|; iTremonti et al 



20041 ). Second. ISearTej (Il97ll ) and lSmithl (|l975h have estal> 



lished many years ago the presence of radial abundance gra- 
dients in the disks of spiral galaxies, with the maximum oxy- 
gen abundance occurring at their centers. Taken together, 
these two facts suggest that the observed oxygen abundance 
should reach its maximal value at the centers of the most 
luminous galaxies. 

In the following, we will be using these notations 
for the line fluxes : R 2 = /[o/i]A3727+A3729 / 'lup, R3 = 



/Ih/3, R.23 = R-2 + 



^[O/I/]A4959+A5007/-^-f//9, R — -f[Oi//]A4363 

R3. With these definitions, the excitation parameter P can 
be expressed as: P = R3/(R2+Rs). The plan of our study 
is as follows. In Section 2 we determine the central oxygen 
abundance for a large sample of nearby spiral galaxies and 
derive the maximum oxygen abundance. The value of the 
oxygen yield is estimated in Section 3. We summarize our 
conclusions in Section 4. 



2 THE MAXIMUM VALUE OF THE OXYGEN 
ABUNDANCE 

2.1 Abundance derivation 

The [OIII]A4363 auroral line is not detected in the majority 
of H 11 regions in spiral galaxies. For a quarter of a century, 
different versions of th e one-dimen s ional empirical calibra- 
tion, first proposed bv IPaeel et al.l l| 19791 ). have b een used 
for abundance deter mination in such H 11 regions. IPilvugin! 
lj200d . l2001al . l2003bh has found that the oxygen abundances 
so determined have a systematic error, depending on the 
excitation parameter defined above. Recently, a relation- 
ship between the observed auroral and nebu lar oxygen line 
fluxes, the ff relation , has been established (|Pilvuginll2005l ; 
IPilvugin et al.|[2006al ). The ff relation allows one to estimate 
the flux in the [OIII]A4363 auroral line, and hence to ap- 
ply the T method to determine abundances in Hn regions 
where that line is not detected, and where only strong oxy- 
gen line intensities are seen. We h ave also derived a ne w 
model-independent ti - t$ relation (|Pilvugin et al.ll2006bh 



t% = 0.72t 3 + 0.26. 



(1) 



We have thus found a way to estimate the electron tem- 
perature in both in the ++ and in the + zones in high- 
metallicity Hn regions where the [OIII]A4363 auroral line is 
not detected. Next, we apply the method described above to 
derive new oxygen abundances in spiral galaxies. Our aim 
is to search for a maximum value of the observed oxygen 
abundance in galaxies. 
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Figure 1. The oxygen abundance as a function of galactocentric 
distance in the disks of the four most luminous spiral galaxies in 
our sample: NGC 1068, NGC 6384, NGC 7331 and IC 342. The 
circles are (0/H)fj abundances in individual Hn regions. The 
lines are linear least-square best fits to those data points. 



2.2 A plateau in the L 
luminosities 



Z diagram at high 



We first define our observation al sample. We use the compi- 
lation bv lPilvugin et all ({2004) of a large sample of published 
spectra of H 11 r egions in nearby spiral galax ie s. Recent mea- 
surem ents from iBresolin et~aT] (1 20041 , 120051 ): ICrockett et~aH 
(2006) have been added to that list. Since the maximum 
value of the observed oxygen abundance is expected to oc- 
cur at the centers of the most luminous galaxies, we first 
derive the radial distribution of the oxygen abundance in 
the disks of the four most luminous spiral galaxies in our 
sample: NGC 1068 with an absolute blue magnitude Mb = 
-22.18, NGC 6384 with M B = -22.22, NGC 7331 with M B 
= -22.20, and IC 342 with M B = -22.27. Th e distances and 
lumin osities of the galaxies are taken from IPilvugin et al.l 
l|2004h . 

Fig. [T] shows the oxygen abundance as a function of 
galactocentric distance in the disks of the four most lumi- 
nous spiral galaxies in our sample. The points are (0/H)ff 
abundances in individual Hll regions. The lines are linear 
least-square best fits to those points. Inspection of Fig. [T] 
shows that the central oxygen abundance in the most lumi- 
nous spiral galaxies can be as large as 12+log(0/H) ~ 8.87. 
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Figure 2. The luminosity — central metallicity diagram. Top 
panel. The open squares show the central (0/H)ff abundances in 
the disks of spiral galaxies plotted against their luminosities. Th e 
filled circles denote irregular galaxies from lPilvugin et al. 
Bottom panel. A blow-up of the high-luminosity part of the L — 
Z diagram in the top panel. 



Are those galaxies indeed the most oxygen-rich ones as ex- 
pected from the L - Z correlation? To clarify this point, we 
next derive the central oxygen abundances of all the spi- 
ral galaxies in our sample and examine their dependence on 
galaxian luminosity. 

The top panel of Fig. [5] shows the familiar L - Z rela- 
tion, where Z denotes here the central oxygen abundance. 
The squares show the central (0/H)g abundances of the 
spiral galaxies in our sample. For the majority of galaxies, 
all H 11 regions with available oxygen emission line measure- 
ments have been used in the abundance gradient analysis. 
For a few galaxies however, some Hn regions have been 
omitted for the following reason. The R23 - O/H relation is 
double- valued, with two distincts parts usually known as the 
"lower" and "upper" branches. The ff relation is valid only 
for Hn regions which are located in the upper branch, with 
12+log(0/H) higher than ~ 8.25. Thus, one has to know 
a priori on which of the two branches the Hn region lies. 
Given that disks of spiral galaxies show radial oxygen abun- 
dance gradients, in the sense that the oxygen abundance is 
higher in the central part of the disk and decreases with 
galactocentric distance, we first consider Hll regions in the 
central part of disks and move outward until the radius R* , 
where the oxygen abundance decreases to 12+log(0/H) ~ 
8.25. The exact value of R* is somewhat uncertain due to the 
scatter in the oxygen abundance values at any fixed radius. 
It should be emphasized that a misuse of ff relation in the 
determination of the oxygen abundance in low-metallicity 
Hll regions beyond R* may produce an artificial change of 



slope at R* in the abundance gradients (|Pilvuginl l2003al ) . 
Therefore, Hll regions with galactocentric distances larger 
than R* , i.e. those with 12+log(0/H) less than 8.25 were 
rejected. 

The filled c ir cles in Fig. [2] show irr egular galaxies from 
iPilvuginl (|2001bh ; |piTvugin et all (|2004T l. Our L - Z diagram 
for irregular galaxies is i n good agreement with the ones 
deriv ed by other authors (jRicher fc McCalllll995l ; iLee et al.l 

2003) . Examination of Fig. [5] shows that the central oxygen 
abundances in galaxies with Mb ~ -20.25 are in the same 
range as those of galaxies with Mb ~ -22.25, i.e. the central 
oxygen abundances in luminous galaxies with Mb < -20.25 
do not show an appeciable correlation with galaxy luminos- 
ity. This flattening of the L - Z relation can be seen more 
clearly in the bottom panel of Fig.[5]which shows a blow-up 
of the high-luminosity end of the relation. The presence of 
an upper envelope in the L - Z relation at 12+log(0/H) ~ 
8.87, is clearly seen. The flattening of the L - Z relation at 
high l u minosities has been n oted before by IPilvugin et al.l 
j2004l ); lTremonti et all (|2004T l. 

What is the meaning of such a plateau? There are two 
known reasons for the existence of the L - Z relation. First, 
it can be caused by a dependence of the efficiency of galac- 
tic winds to get rid of metals on the galaxy's luminosity: 
more luminous and massive galaxies are less efficient in los- 
ing heavy elements by galactic winds. In this case, the L 
- Z relation represents the ability of a given galaxy to re- 
tain the products of its own c hemical evolut ion rather than 
its ability to produce metals (iLarsonl Il974) . It is believed 
that the galactic winds do not play a significant ro le in the 
chem i cal evolution of the largest spiral galaxies (e .g. lGarnettl 
120021 ; IPilvugin et all |2004| ; iTremonti etafl |2004 ). The sec- 
ond reason that has been invoked to explain the L - Z re- 
lation, which we adopt, is that the astration level increases 
and the gas mas s fraction decreases as the luminosity of a 
galaxy increases (jMcGaugh fc de B lok 1997). The existence 
of a plateau at high luminosities in the L - Z relation thus 
implies that, on average, the central parts of large luminous 
spiral galaxies have similar astration levels. 

Another prominent feature of the L - Z relation is the 
rather large scatter in oxygen abundances at a given lumi- 
nosity, Alog(0/H) ~ 0.25. Part of this scatter may be due to 
uncertainties in the oxygen abundances. But another part is 
likely to be real. This scatter may be explained by fluctua- 
tions of the gas mass fraction /1 among galaxies of a given lu- 
minosity. The simple model of chemical evolution of galaxies 
predicts that a decrease of \x by 0.1 results in an increase of 
the oxygen abundance by ~ 0.13 dex, in the range of fi from 
~ 0.50 to ~ 0.05. Then, a scatter Alog(0/H) ~ 0.25 may be 
explained by fluctuations of the gas mass fraction as large as 
A/i ~ 0.2 among galaxies of a given luminosity. The global 
gas mass fractions in the sam ple spiral galaxies have been 
estim ated to be low, /1 < 0.25 (|Garnettll2002l ; IPilvugin" et al.l 

2004) . This suggests that the gas in the centers of the most 
metal-rich galaxies has been almost completely converted 
into stars. Consequently, the observed oxygen abundance in 
the centers of those galaxies represents the maximum attain- 
able value of the oxygen abundance. This provides a natural 
explanation for the constant maximum value of the observed 
central oxygen abundance in the most oxygen-rich galaxies 
(Fig.®. " 

The maximum value of the gas-phase oxygen abun- 
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dance in Hn regions of spiral galaxies is thus 12+log(0/H) 
~ 8.87. Some fract i on of the oxygen is locke d into dust 
gr ains llMever et al. Ill998l ; lEsteban et al. 1998). According 
to lEsteban et aL 1 1998h . the fraction of the dust-phase oxy- 
gen abundance in t he Or ion nebula is about 0.08 dex (but see 
ISimon-Dfaz et al] (|2006l )V Then, the maximum value of the 
gas+dust oxygen abundance in H n regions of spiral galaxies 
is 12+log(0/H) ~ 8.95. 



2.3 Discussion 

The oxygen abundances we have derived here for lumi- 
nous spiral galaxies are significantly lower than those ob- 
tained previou s ly by a number of previous investigators 



tamed previou s ly by a number ot previous inv estigators 
llGarnettl [20021: ITremonti et all 12004 : iMelbourne fc Salzerl 



120021 : lLamareille et al"1l2004r i. This is not surprising. The 



abundances in the papers quoted above have been derived 
using di fferent versions of the early calibratio n. It has been 
argued l|Pilvuginll2001al : iPilvugin et~al1l2004 ) that the oxy- 
gen abundances so determined are signifi cantly overesti- 
mated at the high- met allicity end. Indeed, ITremonti et al.l 
(|2004h have themselves noted that their oxygen abundances 
may have been overestimated by as much as a factor of two. 
On the other hand, the oxygen abundances obtained here 
are not based on any calibration since they are derived via 
the T e method coupled with the ff relation. The ff relation 
is purely empirical in the sense that it relates directly mea- 
sured quantities, without any other assumption. 

The central oxygen abundances of two luminous spi- 
ral galaxies, M 101 (Mb = -21.65) and M 51 (Mb = - 
21.48), have recently been determined from measurements 
of (0/H)t g abundances in a number of Hn regions. They 
have been found to be r espectively 12+log(0/H)T c = 
8.76 dKennicutt et all I200I ). and 12+log(0/H) Tc = 8.72 
|Bresolin et al.ll2004 ). Those data are in agreement with our 
L - Z relation, but are in severe conflict with relations based 
on O/H abundances derived with the one-dimensional em- 
pirical calibrations. 

Comparison with H 11 region photo ionization models has 
led some authors (e.g. Stasihska 2005) to question the ap- 
plicability of the classic T e method to the high-metallicity 
regime. We have suggested (|Pilvugi n 2003b) using the inter- 
stellar oxygen abundance in the solar vicinity, derived with 
very high precision from high-resolution observations of the 
weak interstellar OIA1356 absorption line towards the stars, 
as a "Rosetta stone" to check the reliability of the oxygen 
abundances derived in Hn regions with the T e method. The 
agreement between the value of the oxygen abundance at the 
solar galactocentric distance derived from the T e method 
and that derived from the OIA1356 interstellar absorption 
line towards the stars provides strong support for the appli- 
cability of the classic T c method to solar-metallicity objects. 
We have also examined previously the reliability of oxygen 
abundances derived in high-metallicity H n regions with the 
T e method by analyzing the radial distribut ions of oxygen 
abund ances in the disks of sp i ral ga laxies l|Pilvugin et al.l 
l2006ah . According to lStasiriskal (|2005l l. the derived (0/H) To 
value is very close to the real one as long as the metal- 
licity is low. Discrepancies appear for oxygen abundances 
12+log(0/H) ~ 8.6 - 8.7, and may become ver y large as 
met allicity increases (Fig. la in Stas ihskal [iooH ). The de- 
rived (0/H)t o values are smaller than the real ones, some- 



times by enormous factors. If this is the case, then the ra- 
dial distribution of (0/H)t c abundances should show a bow- 
shaped curve with a maximum value of 12+log(0/H) ~ 8.7 
at some galactocentric distance. However, the derived radial 
distributions of (0/H)t o abundances in the disks of spiral 
galaxies do not show such an appreciable curvature, and 
the (0/H)t c abundances increase more or less monotoni- 
cally w ith decreasing galactocentric distance l|Pilvugin et alj 
l2006al ). Thus, while our results do not rule out the possi- 
ble existence of the Stasihska's effect, they suggest that the 
great majority of Hn regions in galaxies are in a met allic- 
ity range where this effect is not important. Another factor 
that may affect our abundance det erminations is t empera- 
ture fluctuations inside H 11 regions l|Peimbertll 19671 ) . If they 
are important, our abundances would be underestimated. 

The abundances derived here thus depend on the abil- 
ity of the classic T c method to give correct abundances in 
the high-metallicity regime. They will be subject to revision 
if, when metal-rich Hll regions are better understood, it is 
shown that the T c method, because of temperature fluc- 
tuations, strong electron temperature gradients inside Hll 
regions, or any other reason, results in incorrect abundances 
at high metallicities. 

The Sun is one of the widely used reference objects in 
astrophysics. Standard practice is to express the element 
content in a cosmic object via the corresponding value for 
the Sun, i.e. the composition of the Sun is used as stan- 
dard unit. For many years, the recommended solar oxygen 
abundance was 12+log(0/H)Q w 8.9. This high abundance 
was obtained from a one-dimensional hydrostatic model 
of the solar atmosphere. Recently the solar oxygen abun- 
dance has been significantly reduced as a result of a time- 
dependent, three-dimensional hydrodynamical model of the 
solar atmosphere. Taking the average of recent determina - 
tions ( 12+1ok(Q/H)^ =8.70 in lAllende Prieto et al.l |2004l) . 
8.66 in lAsplund et al.1 (|2004h . 8.59 in lMelendej (|2004R the 
solar abundance is now 12+log(0/H)Q = 8.65. Thus, the 
maximum value of the gas+dust oxygen abundance of H 11 
regions in spiral galaxies is higher by a factor of ~ 2 than 
the solar value. 



3 THE OXYGEN YIELD 
3.1 Basic considerations 

The simple model of chemical evolution of galaxies predicts 
that the oxygen abundance of the interstellar matter of a 
galaxy is related to the gas mass fraction fi and the oxygen 
yield Yo by the following formula 



Yo = 



Zq 

Hjt) 



(2) 



In a real situation, the oxygen abundance is also affected by 
the mass exchange between a galaxy and its environment. 
This mass exchange can alter the above relation and mimic 
a variation in the oxygen yield. In that case, the simple 
chemical evolution model is used to estimate the "effective" 
oxyg en yield Y e g (|Edmundslll990l ; IVila-Costas fc Edmunds! 
11992^ 

As noted above, it is believed that galactic winds do 
not play a significant role in the chemical evolution of the 
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largest spiral galaxies. For these, is the oxygen yield derived 
by using Eq. ([2} close to the true oxygen yield? This may 
not be the case for two reasons. First, the simple model is 
based on the instantaneous recycling approximation. Oxy- 
gen is produced and ejected into the interstellar medium by 
massive stars with lifetimes much shorter than the evolution 
time of spiral galaxies. From this point of view, the instan- 
taneous recycling approximation is justified. However, the 
value of Yo depends not only on the amount of oxygen but 
also on the total mass of matter ejected in the interstellar 
medium (see below). In other words, the instantaneous re- 
cycling approximation assumes that the next generation of 
star is formed when all stars from previous generations have 
finished their evolution. This never occurs in a real galaxy. 
As a consequence, the simple model predicts slightly higher 
oxygen abundances than numerical models for the chemical 
evolution of ga laxies with realistic star formation histories 
(Pilvugin 1994). Fortunately, this difference is small and can 
be neglected in the case of oxygen. 

Furthermore, it is well known that the simple model 
predicts many more low-metallicity stars than are observed 
in the solar neighbourhood, the so called "G-dwarf" para- 
dox. Various versions of the infall model, in which an infall of 
gas onto the disk takes place for a long time, have been pro- 
posed to account for the observed m et allicity distribution in 
the so lar neighborhood. ( Tosi 1988a b; Pilvugin fe Edmunds! 

200 ll . among many others) . An in- 



ll996al lbl; [Chiappini et al 



fall model h as also been applied to other spiral galaxies 
|Pagelll 19971 . and references therein). It is thus generally ac- 
cepted that gas infall plays an important role in the chemical 
evolution of disks of spiral galaxies. Therefore, the applica- 
tion of the simple model to large spiral galaxies to estimate 
the true oxygen yields may appear unjustified. It is expected 
that the rate of gas i nfall onto the disk decreases exponen- 
tially with ti me (e.g. Pil vugin & Edmunds 1996a][b|). It has 
been shown (|Pilvugin fc F crrini 1998) that the present-day 
location of a system in the fi - O/H diagram is governed 
by its evolution in the recent past, but is independent of 
its evolution on long timescales. Therefore, the fact that the 
present-day location of spiral galaxies is near the one pre- 
dicted by the simple model is not in conflict with a picture 
in which an infall of gas onto the disk takes place during 
a long time (the latter is necessary to satisfy the observed 
abundance distribution function and the age - metallicity re- 
lation in the solar neighbourhood) since these observational 
data reflect the evolution of the system in the distant past. 
Therefore, one can expect that the application of the simple 
model to large spiral galaxies, Eq. ((2]), provides a more or 
less reasonable estimate of the true oxygen yield Yo- Cer- 
tainly, to find an accurate value of the true oxygen yield 
Yo, an appropriate models of chemical evolution of galaxies 
should be computed. 



3.2 An empirical estimate of the oxygen yield 

From Eq. (J2j| , it is clear that an accurate determination of 
the oxygen yield depends on accurate oxygen abundance and 
gas mass fraction measurements. Our derived (0/H)ff abun- 
dances are more accurate than the (0/H)r 23 abundances 
used in previous studies. However, a precise estimate of the 
gas mass fraction is not a trivial task, and this for several 
reasons. On the one hand, the mass of the stellar component 



of a galaxy is usually estimated by converting the measured 
luminosity to mass via an adopted mass-to-luminosity ra- 
tio. The latter is strongly model-dependent, therefore it is 
difficult to get a reliable estimate of it for individual galax- 
ies. It is widely accepted that the converting factor from the 
near-infrared luminosity is a robust quantity for deriving the 
stellar mass of a galaxy. This quantity shows a low depen- 
dence on the star formation history but depends strongly on 
the initial mass function (for example, on the adopted lower 
stellar mass limit). On the other hand, the mass of molec- 
ular hydrogen can only be estimated by indirect methods. 
The commonly accepted method is the use of the CO line 
flux and a conversion factor X between the flux in the CO 
line and the amount of molecular hydrogen. The conversion 
factor X = N(H2)/I(CO) depends strongly on the physical 
properties of the interstellar medium which are known to 
vary from galaxy to galaxy. The best-estimated values of X 
for a sample of well-studied nearby galaxies s pan the range 
from 0.6 to 10 x 10 20 mol cm" 2 (Kkms -1 )" 1 (|Boselli et all 
2002). Thus, the value of the oxygen yield derived from 
Eq. © can be strongly affected by the uncertainty in the 
gas mass fraction determination. 

The use of the maximum value of the oxygen abundance 
derived here allows us to overcome the above problem in the 
following way. We compare the derived maximum value of 
the oxygen abundance in galaxies with oxygen abundances 
predicted by the simple models with different oxygen yields 
for jj,—0. Fig. [3] shows the oxygen abundance as a function 
of the gas mass fraction predicted by the simple model with 
Y o =0.0030 (solid line) and with Y o =0.0035 (dashed line). 
Since the oxygen abundances are expressed in units of num- 
ber of oxygen atoms relative to hydrogen, while Zo in Eq. @ 
have units of mass fractio n, we adopt th e following conver- 
sion equation for oxygen (|Garnettll2002T ) 



Zo = 12 



(3) 



The simple model breaks down as the gas mass fraction 
approaches zero because the term ln(l//i) blows up. There- 
fore the oxygen abundance predicted by the simple model 
at fi—0 should be estimated by extrapolation of the model 
predictions to fi=0. The extrapolations of the simple models 
for two values of the oxygen yield are shown in Fig. [3] by the 
dotted lines. The fact that the central oxygen abundances 
in galaxies are derived here, not as oxygen abundances of 
Hn regions in the very central parts of galaxies, but also as 
extrapolations of linear fits, justifies the above method. 

The maximum value of the observed gas-phase oxygen 
abundance in Hn regions of spiral galaxies is shown by the 
open square in Fig. [3] The maximum value of the gas+dust 
oxygen abundance is shown by the open circle. Examination 
of Fig.|3]shows that the maximum value of the observed gas- 
phase oxygen abundance in Hn regions of spiral galaxies 
corresponds to the simple model with Yo ~ 0.0030, and the 
maximum value of the gas+dast oxygen abundance in spiral 
galaxies corresponds to the simple model with Yo ~ 0.0035. 

Thus, we can conclude that the value of the oxygen 
yield is about 0.0035, depending on the fraction of oxy- 
gen incorporated into dust grains. The value of the oxygen 
yield derived here from the gas-phase oxygen abundance 
(Yo ~ 0.0030) is close to tha t obta ined recently for spi- 
ral galaxies by IPilvugin et ail \2004 ) (Y Q « 0.0027) and 
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Figure 3. The oxygen abundance as a function of the gas mass 
fraction predicted by the simple models of chemical evolution of 
galaxies, with an oxygen yield Yo=0.0030 (solid line) and with 
Yo=0.0035 (dashed line). The extrapolations to fi=0 for these 
models are shown by the dotted lines. The open square is the 
maximum value of the observed gas-phase oxygen abundance in 
H II regions of spiral galaxies. The open circle is the maximum 
value of the gas+dust oxygen abundance. 



bv lBresolin et alJ l|2004 ) (Yo « 0.0032) but is s i gnific antly 
lower than the oxygen yield obtained bv lCarnettl ( 2002I ) (Yo 
» 0.010). 

3.3 The stellar oxygen yield 

Can the derived value of Yo be reproduced by the ex- 
isting stellar evolution models? The similarity between 
the present oxygen abundances in the interstellar mat- 
ter in the solar vicinity and that in the Sun which was 
set up some 4.5 Gyr ago, suggests that spiral galax- 
ies have had a high oxygen abundance during most of 
their evolution. We therefore use estimates of the yields 
in the high-metallicity regime. We have compi led the 

(IMaederl 



19921; Wooslev & Weaver 19951; 


Langer & HenkeT 19951; 


Nomoto et al.ll997; Portinari et al. 


1998;lMevnet & Maederl 


2002: Hirschi et all 


2005 


: iKobavashi et all 120061). Follow- 


ing iTimmes et al.l ( 


1995 


), we have used models A from 
) for stars in the 11-25 Mq mass 


Wooslev & Weaver! ( 


1995 



range, and models B for masses equal to 30, 3 5, and 40 Mq. 
The yields computed by iMevnet fc Maederl (l2002j ) for th e 
metallicity Z = 0.020 are given in lChiappini et al.l (|2003al ) . 
The mass Mo of freshly produced oxygen ejected by a star in 
the interstellar medium, calculated by different authors, is 
shown as a function of the initial stellar mass in the top 
panel of Fig. [4] That figure shows that there exists sig- 
nificant differences between the various predictions of Mo 
at a given initial stellar mass. The differences are espe- 
cially large for stellar masses larger than ~ 25 Mq. They 
are caused by differences in the input physics (stellar wind, 
mixing, etc) used by the various authors. Close examina- 
tion of the top panel of Fig. U reveals that there are two 
main gr o ups o f points: a group of high Mq values from 
IMaederl l| 19921) ; I Wooslev fc Weaver! (| 19951 ); iNomoto et ail 



jl997l ); lKobavashi et all (2006), based on stellar models with 
small mass loss r ates from stella r winds; and a grou p of low 
Mp values from IMaederl l|l992l ); lLanger fc Henkell (|l995h ; 
iPortinari et alJ (| 19981 ). based on stellar models with large 
mass loss rates from stellar winds. We will consider both 
the high and low values of Mo- They can be considered as 
delimiting the range of oxygen production by stars. 

It worth noting that the oxygen production is computed 
for single stars. A large fraction of stars are members of bi- 
nary systems and evolve with mass exchange between com- 
ponents. Therefore one may expect that the stellar evolution 
models with strong stellar winds give values of the oxygen 
yield that are closer to the true value. 

The cumulative mass Qo of freshly manufactured oxy- 
gen ejected by a single stellar population at time r after 
its formation is calculated as the sum of contributions from 
each star down to a stellar mass Mg, corresponding to a 
lifetime r. It is given by the following expression 



Qo{M* s ,M v ) 



E 

Mi 



Mq if A Ms 



(4) 



where ip AMs is the number of stars within the mass interval 
AMs as determined from the initial mass function if, and 
Mo is the ejected mass of freshly manufactured oxygen from 
a star with mass Ms. Since stars with masses lower than 10 
Mq do not manufacture an appreciable amount of oxygen, 
the lower mass limit Mg in Eq. Q has been set to 10 Mq. 
Difficulties arise in the computation of Qo because of the 
limited number of stellar masses for which models exist. The 
gaps between the masses for which models have been calcu- 
lated produce artificia l jumps in the resulting Qo(MtjMs) 
l|Pilvuginl 1 19921 . 1 19941 ). To overcome these difficulties, we 
have adopted the following approach. The dependence of the 
ejected mass of freshly manufactured oxygen Mo on stellar 
mass Ms is approximated by a polynomial of degree 3. We 
have adopted the approximation 



log Mo = -23.29 + 38.84 log M s 

-21.02 (log M s ) 2 + 3.84 (logM s ) 3 



(5) 



for the case of models with no or low stellar wind mass loss 
rates. This relation is shown by the solid line in the top panel 
of Fig. U As for the case of models with high stellar wind 
mass loss rates, we adopt the approximation 



log Mo = -24.04 + 41.69 log M s 

-23.66 (log M s ) 2 + 4.46 (logM s ) 3 



(6) 



This relation is shown by the dashed line in the top panel 
of Fig. [4] Model points with large deviations were not used 
in the derivation of these two approximations. 

The dependence of the mass M re t of the matter returned 
by a star to the interstellar medium on the initial mass of 
the star is shown in the bottom panel of Fig. U Inspection 
of the figure shows that that there is agreement between the 
Mr C t values from different authors, so that all the data can 
be approximated by the following expression 



log M rct = - . 344 + 1 . 666 log M s 

-0.500 (logMs) 2 + 0.126 (logM s ) 3 



(7) 



The cumulative mass Q TCt returned by a single stellar pop- 
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Figure 4. Top panel. The mass of freshly manufactured oxy- 
gen ejected by a star in the interstellar mediu m as a function 
of the initial mass of the star. The results from [Maederl lll992h 
for models with low stellar wind mass loss rates are shown 
by open squares, and that for models with high stellar wind 
mass loss rates are shown by filled sq uares. Results from other 
investigators have also been plott ed: lLanger fc Henkell dl995h 
(filled circles), iNomoto et alj lll997l) fplus signs) , I Portinari et al.l 
U998) (filled triangles). iMevnet fc Maederl J2002T) (open rhom- 
bus), Hirschi et alj ll2005l) (filled rhombus), and Kobavash i et all 
(2006) (o pen triangles). The solid line is a fit to the model 
results of iMaederl lll992l) (models w i th low stellar wind mass 
loss rates l.lWooslev fc Weaverl lll995l) ,l Nomoto et al. 

I <1997l) . and 

IKobavashi et al l (1200^ 7 The dashed line is a fit to the model re- 
sults of iMaederl (Il992l) ( models wit h high stellar wind ma ss loss 
rates). lLanger fe Henkel dl995h and IPortinari et al.l lll998l) . Bot- 
tom panel. The mass returned by a star to the interstellar medium 
as a function of the initial mass of the star. The data for low- 
and in termediate-mass stars from Ivan den H oek fc Groenewegenl 
(1997) are shown by crosses, those from iMarigol ll200ll) by open 
circles. The sources for the models and the symbols for massive 
stars are the same as in the top panel. The solid line is a fit to all 
the data. 



ulation is given by the following expression 

Qrct (Mg , Mv) = ^ M rct ipAM s . 



(8) 



To compute Qo and Q ro t , we have to p recise the initia l 
mass function (IMF) . We first consider the ISalpeteJ (|l955l ) 
IMF 
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Figure 5. Top panel. The lower mass limit Ml of the Initial Mass 
Function (IMF) calculated so that, for a given upper mass limit 
Mtj, the oxygen yield Yo is equal to 0.0035. The circles corre- 
spond to Kroupa's IMF and oxygen yields given by stellar models 
with high stellar wind mass loss rates. The squares correspond to 
a Salpeter IMF and oxygen yields given by stellar models with 
high stellar wind mass loss rates. The triangles correspond to a 
Salpeter IMF and oxygen yields given by stellar models with low 
stellar wind mass loss rates. Masses are in units of solar masses. 
Bottom panel. The fraction f of the total stellar mass in stars with 
masses above IMq for the three combinations of IMF parameters 
and oxygen yields in the top panel. 



The coefficient c is defined by normalizing the IMF over the 
whole mass interval 



Mu 



M s 0(M s )dMs 



(10) 



<MM £ 



cMo 



(9) 



where Ml and Mu are the lower and upper mass limits. 
These mass limits are not known and are in fact free pa- 
rameters. For example, Romano et alj (|2005l ) adopt Mu = 
100M Q while IKobavashi et alj i|2006t ) choose Mu = 5OM . 
Sometimes, the fraction ( of the total stellar mass locked 
up in stars with mass es above 1M© is fixed, instead of Ml- 
IPortinari et al.1 (|l998l ) have found that a "good range" for £ 
is between 0.3 and 0.4. 

We next check whether there exist values of Ml and Mu 
so that Yg od = Yo, where the model oxygen yield Yg od is, 
by definition, the ratio of the mass of newly manufactured 
oxygen ejected by a single stellar population to the mass 
locked up in long-lived stars or remnants of that single stellar 
population 
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y-mod 
o 



MM') 



•Qret(Af*) 



(11) 



Since stars with masses less t han 0.95 Mp have lifetimes 
longer than the Hubble time (|Schaller et al.l [l992h . we set 



M* 



0.95 M in Eqs. (JUTTJ . Then, by fixing Yg° 



0.0035 



Cl 


Ma 1 - 3 


for 


M s < 0.5 M 




C2 


Mg 2 ' 2 


for 


0.5 M Q < M s < 1M© 


(12) 


C2 


M s " 2 - 7 


for 


M s > 1 M© 





and Mu, we can find Ml. The resulting dependence of Ml 
on Mu is shown in the upper panel of Fig[S] The triangles 
correspond to oxygen productions of stellar models with low 
stellar wind mass loss rates (Eq. (JSJl), while the squares cor- 
respond to oxygen productions of stellar models with high 
stellar wind mass loss rates (Eq. l[6)l). The lower panel of 
FigE] shows by the same symbols the fraction C, of the total 
stellar mass in stars with masses above 1M©, for the same 
combinations of IMF parameters and oxygen productions as 
in the upper panel. 

A number of multi-component power-la w stellar IMFs 
have been proposed as well. The one from IKroupa et al.l 
i|l993h is often used in the construction of chemical evo- 
lution models of galaxies. IRomano" et al.1 (120051 s ) have found 
that the Kroupa IMF fits better several observed properties 
of the solar vicinity as compared t o the Salpeter IMF. The 
three-component power-law IMF of lKroupa et al.l (1 19931 ) can 
be parameterized as follows 



0(M S ) 



The coefficients Ci and C2 are defined by normalizing the 
IMF over the whole mass interval (Eq. (|10[) ). and by the 
condition ci Mg 13 = c 2 Mg 2 ' 2 for M s = 0.5M©. 

The dependence of Ml on Mu for the Kroupa IMF and 
for oxygen productions from models with high stellar wind 
mass loss rates (Eq. is shown by circles in the top panel 
of Fig[5] The bottom panel shows also by circles the fraction 
£ of the total stellar mass in stars with masses above 1M© , 
for the same IMF parameters and oxygen productions as in 
the top panel. The value Yo = 0.0035 cannot be reproduced 
with a Kroupa IMF and oxygen productions from models 
with low stellar wind mass loss rates. 

In summary, we have found that in the case of a Salpeter 
IMF, the observed yield, Yo = 0.0035, can be satisfactorily 
reproduced by existing stellar evolution models with both 
low and high stellar wind mass loss rates. However, in the 
case of a Kroupa IMF, the observed yield can only be re- 
produced by stellar evolution models with high stellar wind 
mass loss rates. If the Kroupa IMF is more realistic than 
the Salpeter IMF, then our results suggest that stellar evo- 
lution models with low stellar wind mass loss rates predict 
too large oxygen yields. In the case of stellar models with 
high stellar wind mass loss rates, we obtain £ ~ 0.25 for a 
Salpeter IMF and ( ~ 0.40 for a Kroupa IMF. These val- 
ues are close to the " good range", 0.3 < C < 0.4, found by 
iPortinari etaH l|l99Sl) . 

The above conclus ion is at odds with the finding of 
IChiappini etU] (l2003bh that the oxyg e n yie lds in massive 
stars computed by IWooslev fc Weaver! (1 19951 ) without tak- 
ing into account mass loss reproduce well abundance obser- 
vations in t he spiral galaxy M 101. The origin of the di- 
vergences m IChiappini etal] l|2003bl ) conclusions and ours 
probably comes from the fact that those authors used an 
early higher value of the central oxygen abundance of M 101, 
12+log(0/H) ~ 9.2, while recent measurements of (0/H)t o 



abundances in a number of Hn regions in M 101 have 
given a central oxygen abundance 12+log(0/H)T c = 8.76 
(|Kennicutt et alj|2003h . 



4 CONCLUSIONS 

We search here for the maximum oxygen abundance in spiral 
galaxies. It is expected that this maximum value occurs at 
the centers of the most luminous galaxies. The luminosity - 
central metallicity diagram for spiral galaxies is constructed. 
The central oxygen abundance in a galaxy is derived from a 
linear least-square best fit to the abundances in individual 
Hn regions at various galactocentric distances. The oxygen 
abundance in an H n region is derived using the T e method, 
with the flux in the [OIII] A4363 auroral line determined from 
the ff relation when not available from observations. The 
electron temper ature t2 is estimated from a newly derived 
t2 - t3 relation l|Pilvugin et al.ll2006bl ). 

We found that there exists a plateau in the luminosity 
- central metallicity diagram at high luminosities (-22.3 < 
Mb < -20.3). This provides strong evidence that the oxygen 
abundance in the centers of the most metal-rich luminous 
spiral galaxies reaches the maximum attainable value of oxy- 
gen abundance. The maximum value of the gas-phase oxygen 
abundance in Hn regions of spiral galaxies is 12+log(0/H) 
~ 8.87. Because some fraction of the oxygen (about 0.08 
dex) is expected to be locked into dust grains, the maxi- 
mum value of the true gas+dust oxygen abundance in Hn 
regions of spiral galaxies is 12+log(0/H) ~ 8.95. This value 
is a factor of ~ 2 higher than the recently estimated solar 
value. 

Based on our derived maximum value of the oxygen 
abundance in spiral galaxies, we have estimated the oxygen 
yield. We have found that the oxygen yield is around 0.0035, 
depending on the fraction of oxygen incorporated into dust 
grains. 
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